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Some fifteen years ago, two independent accounts of the concept of planar chirality were proposed in [TJ Ch. 8] 
and [2J, seeking to generalize the notion of conventional, i.e., 3-D chirality [31 0] to 2-D space. In fact, in [HIS], 
2-D chirality arose as one particular class from a reconsideration of the notions of geometrical dissymmetry [3j [6] 
and polarization of electromagnetic (EM) waves in a multi-dimensional context. It lead to the suggestion that 
planar chiral structures allow for the engineering of a novel class of thin-film artificial composite bianisotropic 
metamaterials exhibiting optical activity, which was the primary motivation for that study. 

The concept of 2-D chirality is easy to grasp: when viewed from "above" its plane, a planar 2-D handed 
structure viewed looks like the optical mirror image of the same structure viewed from below the plane. This is in 
marked contrast with 3-D chirality, whose handedness is an intrinsic property of the structure independent of the 
direction of observation (wave incidence): for example, the screw sense of a given helix, which is the archetype 3-D 
chiral structure, is the same when viewed from either end. This demonstrates that spatial dispersion is inherent to 
planar chirality. 

Features of planar chiral structures that are of major practical importance include their ability to be fabricated 
and deposited as layered thin films, their application as magnetic thin films whose properties arise from magneto- 
electric coupling, and their potential for generating unusually large optical activity (i.e., relatively large rotary 
power for the plane of polarization of an impinging wave) compared to that for 3-D chiral structures, as was 
numerically predicted early on [7] and experimentally verified, cf., e.g., [5]-|14|. This rotary power results from the 
fact that, like 3-D chiral geometries, 2-D chirality gives rise to a different sensitivity of incident waves exhibiting 
left- vs. right-hand circular eigenpolarizations (circular birefringence). The rotary power, and its extension to 
diffracted beams, has formed the basis for attracting wider interest for several decades, cf., e.g., [5]-|19j. 

To date, the numerical analyses of planar chiral structures have mainly been performed using full- wave 
numerical methods [TJ |20l [21] , using a formulation of the method of moments in terms of subdomain basis function 
expansions for straight conducting wire segments. In these studies, the main focus was on the calculation of optical 
rotation and ellipticity in the scattered field and on microscopic multi-polarizabilities of planar chiral structures. 
Subdomain formulations provide flexibility for treating arbitrary types and shapes of planar chiral geometries. In 
this paper, instead, we use entire-domain basis functions that are computationally more efficient, more accurate, 
and - although less versatile - better suited to the particular class of planar chiral geometry of interest here, 
viz., asymmetrically split rings. Furthermore, the closed-form expressions for the induced current afforded by an 
entire-domain formulation offer insight into the characteristic EM features of planar chirality, and the role of the 
geometry parameters and impedance loading herein. 

This paper is organized as follows. In Section [21 after giving a motivation for the method, we develop 
the Fourier expansion for the case of a single active multi-split loop with localized voltage sources across one 
or more of its gaps and for a passive loop illuminated by an externally incident plane wave. We show how 
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the electromagnetic characteristics of planar chirality transpire from the expressions for the azimuthal current 
distribution in enantiomeric forms of the loop. The results are valid for abritrary frequencies. In Section [31 we 
use a quasi-static approximation (for use with effective- medium theory for the homogenized medium). Again, 
characteristics of planar chirality are determined, now for the polarizabilities of the loop. In Section (4[ we turn 
attention to diffraction gratings of multi-split loops and determine their reflection and transmission characteristics as 
a function of frequency, using the periodic method of moments. The effect of planar chirality on these characteristics 
is demonstrated. Finally, conclusions of this theoretical study are given Section [5j 

2. Current distribution in multi-split planar chiral ring 

2.1. Motivation 

Previous studies of planar chiral geometries have focused on straight or curved open-ended wire-type geometries, 
such as generalized swastikas (gammadions) [7], [10]-[I7], [21], planar chiral hooks (L-or J-shaped wires, or indeed 
most letters of the alphabet) [TJ [31 [TB] , scalene triangles [331 [22], spirals [MUM], etc. These works follow a "bottom- 
up" approach, in which a planar chiral geometry is synthesized via assembly of concatenated 1-D wire segments. 
However, planar chirality is inherently a 2-D notion, hence it is natural to consider a canonical 2-D geometry formed 
by a conducting loop with unequally spaced incisions (gaps) along its circumference, producing arcs of unequal 
lengths. This "top-down" approach has the advantage of allowing for EM properties to be calculated efficiently 
based on the theory of loop antennas, by expanding the loop current in a Fourier series of current eigenmodes of 
the corresponding closed loop. As will be shown, this analytical method enables us to demonstrate in mathematical 
terms the essence of planar chirality. 

The connection between triple-split rings and triangles [22] is obvious: the locations of three splits in an 
otherwise closed ring define the vertices of a triangle inscribed or circumscribed by the split ring. If their azimuthal 
locations define a scalene triangle then the split ring exhibits planar chirality; otherwise (i.e., if the associated 
triangle is isosceles or equilateral) the split ring exhibits one or three planes of planar symmetry, respectively. The 
top-down approach for studying the chiraity of split loops has a number of advantages over triangles, viz., allowing 
for an elegant analytic framework for the current distribution and EM properties, as well as allowing for effortless 
extension to multiple splits that define inscribed or circumscribed planar polygons. 

Knowledge of the distribution (mode) of the impressed or induced current in the split ring is fundamental in the 
calculation of electric and magnetic dipole and higher-order multipole moments and polarizabilities of planar chiral 
particles. The latter, in turn, form the basis for macroscopic characterization of the permittivity and permeability 
of the homogenized composite material as an effective medium and for its radiation, scattering, and reflection 
properties. By considering each gap in turn as either a generator or a load, a multi-split loop can be conceived as 
a multiply-driven and multiply-loaded loop. 
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2.2. Active ring driven by localized voltage source 

A ring driven at one or more of its gaps by an impressed voltage source allows for constructing an active array whose 
spatial and spectral radiation characteristics are specified. For example, active surfaces combined with double-split 
ring arrays have been proposed for a semi-classical implementation of a spaser [26j . It is conceivable that further 
integration and miniaturization of such arrays may be achieved by using an active element within the gaps of the 
split ring. 

2.2.1. Closed ring Consider an unloaded round- wire loop of radius a located in free space with intrinsic impedance 



Vo = v Mo / e o — 1207T ohm. The wire radius ro is chosen sufficiently small compared to a (i.e., with King index 
k = 21n(27rro/o) larger than 10) so that the loop can be considered in its thin- wire approximation, i.e., the radial 
current component is negligibly small compared to the azimuthal current 1(4'). This current, when produced by a 
driving emf or voltage source V located at azimuth 4 — 4>o, i.e., V5(4 — (f>o), can be expanded in a Fourier series 
as [27], PHI Ch. 9] 

1(4) = Vy(cf>), (1) 

where 

is the local admittance of the loop as measured at 4- In particular, the input admittance at 4o = is y(0) = Y. A 
time-harmonic dependence exp(jcji) for V and I is assumed and suppressed throughout. The Fourier coefficients 
of I (4>) are given by 
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are Lommel-Weber and Bessel functions, respectively, of the first kind and order n. Explicit analytical 
approximations for A n valid for ka < 1.3 can be found in |29j . 

Figure Q] shows typical values of the first twenty coefficients l/a n for a thin-wire loop at ka — 1.45 calculated 
from the previous expressions based on numerical integration. 

1/a n for ka=1.45 




Figure 1. Real (blue) and imaginary (red) components of l/a n (n = 0, . . . , 19) for a closed loop with loop 
radius a = 20.5 mm, wire radius ro = 0.82 mm and electrical size ka = 1.45. 



2.2.2. Multi-split ring Now consider a loop that exhibits N gaps (splits) of respective permittivities and 
widths d(f>k that are centered at arbitrary azimuthal locations <fi — tpk (fc = 1, . . . , N) (Figure^. For simultaneous 
excitation of the loop at all its gaps, the total azimuthal current distribution is found by considering each (fcth) gap 
in succession as a discrete voltage source with the N — 1 remaining gaps presenting local capacitive load impedances 
Zi5((f> — (pi) — S (pi /(ju) efc7rr"o) in a lumped-circuit representation of the loop for this fcth partial excitation, with 
i = 1, . . . , fc — 1, k + 1, . . . , N. Inserting lumped-element impedances or more complicatd impedance circuits across 
the gaps allows for devising more general complex values of Zi with user-defined zeroes and poles in the complex 
plane as a function of frequency, including active impedances and circuits (Ke(Z) < 0) carrying their own voltage 
or current source. 

For the current Ik generated by a delta-distributed voltage source across the fcth gap, 

N 

h{<f) = v kV {(p - 4> h ) - zM<t>i)y{4> - (9) 
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(a) 




Figure 2. (a) Three-split ring with loop radius a, wire radius ro, azimuthal gap locations <f>k, and gap 
widths 5<j>k (k = 1,2,3). (b) External illumination of three-split ring by plane wave. The Cartesian field 
components E' cos if) and E 1 sin 1/) are the projections of the incident electric field vector E_' onto the oxy- 
plane of the ring and on the plane (j> = 4>o, respectively. The angle of incident polarization ip is measured 
with reference to the former component. The elevation angle of incident propagation 9 is measured between 
the normal oz to the plane of the loop and the line connecting the centre of the loop o with the distant 
transmitter, along the direction of incidence j_ fci = fe l /||fe l ||. The positive direction of the azimuthal angle 
of incident propagation c/>o is from the positive oi-axis ((f) = 0) to the projection of the connecting line 
onto the eery-plane. 

Upon summation of the Ik{4>i) with respect to all N — 1 gaps, i.e., summing 

N 

hi<Pi) = v ky (& -<f> k )+Yl Z M<t>M<i>i ~ <t>i) (10) 
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with respect to i for fixed k, this leads to 
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Repeating the calculation for each /cth gap source in succession, this yields a system of N equations in the N 
unknowns Ik{4>i)- Its solution can be substituted back into ([9]) to yield the fcth partial current distributions ifc(0). 
Finally, the total current associated with the simultaneously excitation of all N gaps is then 

N 

i(<t>)=j2u<t>)- (12) 
fc=i 

2.3. Manifestation of planar chirality in azimuthal distribution of current 

The electromagnetic manifestation of planar chiral properties can most readily be demonstrated from the mirror 
symmetries of I{<f>) in a two-split ring with unequal gap widths (Z\ ^ Z 2 ) or in a three-split ring with equal gaps 
(Z\ = Zi = Z%; a fortiori for three unequal gaps or load impedances). We refer to both cases as electromagnetic 
(or load-induced) chirality and geometry-induced planar chirality, respectively. Since both types are determined, 
respectively, by the characteristics of gaps and loop segments totalling 2n, they are in a sense complementary 
(dual) for unloaded gaps. 

For the two-split ring, the partial current distributions are 

4(0) = V k y(<p - fe ) - Z s I k [fl>M<t> ~ &)> (13) 

where {j, k} = {1, 2} with j ^ k. Together with the partial gap currents Ik(<t>j) = VkYjy((f>j — <pk)/[Yj + 2/(0)] an d 
1/Zj =Yj, we obtain 

7(0) = V 1 y( ( f > -<l> 1 ) + V2y{4>-4>2) 
2/(^2 - <pi)y(<p~ 2 ) 
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Y 2 + 2/(0) 
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The current distribution of the corresponding mirror ring is obtained by a reflection in its plane with respect to 
an orthogonal plane at the bisecting azimuth (f>o = (0i + 02)/2, i.e., replacing Vi,Zi by V2,^2 or, equivalcntly, 
replacing 0i, 2 by -0i, -0 2 , yielding 

/mirr(0) = V lV {(j) + 0i) + F 2 J/(0 + 2 ) 



Y 2 + 2/(0) 

- 0i)2/(0 
Yi + 2/(0) 
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y(fa - <f>i)y((t> - fa) 



V 2 



Yi + j/(0) 



T/ y(fa -fa)y(<f>-fa) n „, 

~ Vl y^TW) ' (16) 

The expressions (fT5|) and (fT6|) are equivalent because </>i — fa = — (fa — fa). When now replacing both Vi,Zi 
by V2,Z2 and fa, fa by — (pi,— fa, it is verified that I m i T r(fa transforms back to the original current distribution 
I(fa, thus demonstrating planar chirality. [Note that y(4> + A^/2) ^ y(4> — A0/2) (except at <fi — and 7r) where 
A0 = (fa — fa)/2, because y(-) is an even function.] A canonical case for the two-split ring that interacts with 
both TE (or s) and TM (or p) polarizations arises for the special case where the splits are located at right angles 
(Ac/) = ±7r/2). In the latter case, y(+n/2) = y(-ir/2). 
For the three-gap ring (Figure [2]) , the total current is 

I(fa)=h(fa)+h(4>)+h(4>) (17) 



where 



h(<j>) = V iy (<t> - fa) ~ Z 2 h(fa)y(<p - 2 ) 

-ZshfaMf-fa), (18) 



with 



VIM J \ 1 _ yW2-V3f 

h(fa) - V\ v ^ (n) yUs-toMto-M (19) 

y 2 (0) 

while expressions for i2(</>) and I-i(fa) follow from (|18|) — (f20|) by cyclic permutation of indices. Defining fa — as the 
reference azimuthal orientation, it is found that for a planar achiral structure (fa — — fa) with symmetric (V" 2 = V3) 
or anti-symmetric (V 2 = — V3) excitation, I (fa is a function of y(fa) and y(fa- fa)±y(fa- fa) = y(fa-fa)±y((j>+fa) 
only, viz., 

I(fa = V lV (fa + V 2 [y(cp - fa) ± y(cp + fa)} 

T , y(fa) 1 y+2/(o) r , , , s . /j, 1 j. m 
~ F ' v± ^ — j/W>- <fe) + 2/(0 + ^2) 

T/ y(202j 1 y+y(o) , , , ,,,,,,,,, 
y + 2/(0 i_ yw± 

y 2 (0) 

in which the upper and lower signs refer to the cases V% = V3 and V 2 = — V3, respectively. By contrast, for a planar 
chiral structure, the dependence on y(<f> — fa) and y(4> — fa) does not arise from their sum or difference. 
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2.4- Passive split ring illuminated by incident plane wave 

Passive split rings are the building blocks of metamaterials, metafilms and diffraction gratings, including high-Q 
and/or planar chiral spatial polarizers or filters, as well as electromagnetic bandgap (EBG) materials and disordered 
artificial media. Analysis of a unit-cell single inclusion is therefore of fundamental interest. 

For an incident plane wave with polarization angle tp measured in the transverse plane and with elevation 
angle of incidence 9 (Figure [2Kb)), the normalized admittance of an unloaded closed loop is Ch. 10] 



(<£) = — Y] — exp(-jn0) 



7r?7o a r , 

' n— — oo 



(22) 



in which a n is defined by §3§ and the Fourier expansion coefficients for the incident electric field are given by 

_ j"exp(-jn0 o ) 
h - 2 

x {sin ip cos 9 [J n +i (ka sin 9) + J„_i(fcasin#)] 

+j cos ip [J n+ i(kasin9) — J„_i(fcasin#)]} . (23) 

In particular, for grazing incidence of a TE plane wave whose plane of incidence coincides with the plane of the 
loop itself (9 = 90°, ip = 0°), these coefficients read 



f n =P nexp(-jn^o) 



J n (ka) 
ka 



(24) 



whereas for normal incidence perpendicular to the plane of the loop (9 = 0°, = 0°), the coefficients are 

0. (25) 



f n f exp(-j0 o ) 

JO — U, Ji — , Jn>2 



For an externally illuminated one-split loop containing a single admittance load Yl = 1/^l across its gap at 
= </>l, the current expansion is now 

2/(0 - 0l)" 



i{4>) = v 



u{<P) 



Y + Y h 



(26) 



where V = — C a E\ is the emf induced by the external electric field E\ in the corresponding unloaded ring with 
circumference C a = 2na. This current is to be compared to 

y(M y{<t>- <t>i) 



i(4>) = v 



(27) 



Y + Y L 

for a loop locally driven at <f> — by a gap voltage source V in the absence of external illumination. By extension 
of (f2"6")l , for N loads at N splits located at tfij, the induced current is 



JY 



(28) 



With 



for i = l,. 



N 

I(<f>i) = Vu(^) - Z i I i ( t>M ( t>i - 4>i) (29) 
j'=i 

, AT, the solution of the TV equations (|2^)) in the N variables I(4>i) can be expressed in matrix form as 

-l 



V [I + A, 



(30) 
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where 7 is the unit N x N dyadic and 

I L = [7(0 1 ).../(0 A r)] T (31) 

Ul =[«(&)... tt(^r)] T (32) 

A h = [y(^i-<l>j) Zj]ij=i,...,N- (33) 

The thus obtained current components [7 L ]j = I(4>j) are then substituted back into (|28p to yield /(</>). A locally 
driven ring follows as a special case, by replacing u(4>) with y(4>) and omitting the self-terms j = k by defining 
Zj = Zk per chosen gap j. 

In particular, for a one-split ring, 

7(0) = Vu(<f>) ~ Z L I(fa)y(<t> - <£ L ) (34) 

or, with /(<^l) = VY L u(fa)/{Y L + y(0)} and 1/Z, = Yj, 

7(0)=Vu(0)-V — — — 7- — . (35 
*L + 2/(0) 

For a two-split ring, 

I{4>) = Vu{4>) - - <£i) - Z 2 I(<f> 2 )y(<t> - <j> 2 ) 

(36) 

with 

7(0i) = Vx 

[1 + ^2j/(0)] u((j)i) - Z 2 y((f>l - (f>2)u{(/)2) 

[1 + Z iy {0)} [1 + Z 2 y{0)} - Z iy {cf> 2 - 0i)Z 2 y(0! - fa) 

(37) 

= v x 

[1 + Z iy (0)} u(<h) - Z iy (fo - 4>i)u(h) 



[1 + Z iy (0)] [1 + Z 2 y(0)] - Z lV {<h - <t>i)Z 2 y{<t>i - <f> 2 ) ' 

(38) 

Figure [3] shows 7(0) for a three-split ring with a = 20.5 mm, tq = 0.82 mm and characterized by either 
(<t>i,<h,<fo) = (45°, 90°, 180°) or its planar chiral enantiomer (<f>i,fa,<f> 3 ) = (315°, 270°, 180°), with Sfa = 6(j> 2 = 
(503 = 3° in both cases. The mirror dissymmetry of the structure with respect to the ox-axis causes a mirror 
symmetry of the azimuthal current distribution with respect to <f> = 180°. Note that, at this relatively high 
frequency, deviations from the quasi-static result Re[7(<^)] = — jlm[7(0)] occur, although these deviations are equal 
in magnitude for both mirror image rings. 

Alternatively, for given fixed gap locations, manifestation of planar chirality also occurs through a choice of 
different gap widths but whose values progress monotonically in either clockwise or counterclockwise direction. The 
dissymmetry of the azimuthal current distribution exists for any direction of polarization and incidence, i.e., for 
any value of 9 and ip, including normal and grazing incidence as extreme cases, although it is more pronounced in 
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Figure 3. Real (blue solid curves) and imaginary (green dashed curves) parts of the current /(</>) induced 
in a passive three-split loop with ka = 2, a — 20.5 mm, ro = 0.82 mm, calculated from the first twenty 
coefficients l/a n (n = 0, . . . , 19) for the closed loop and the coefficients /„ associated with a plane wave 
at normal incidence and polarized along the s-direction, i.e., (8, cj>o,ip) = (180°, —90°, 0°): (a) for a planar 
chiral split ring with (c/>i, c/>2, </>3) = (45°, 90°, 180°), 5</>i ~ S(/>2 = <5c/>3 = 3°; (b) for its planar enantiomer 
form with (cfn,^,^) = (315°, 270°, 180°), 50i = Sfo = Sfo = 3°. 



the case of oblique (glazing) incidence as presented here. Figured] illustrates this with the case of a three-split loop 
at ka = 1.5 for (0 = 30°, 4> = 0°,ip = 0°) with (^,^2,^3) = (1°, 10°, 19°) vs. (^1,^2,^3) = (19°, 10°, 1°) 
with identical gap centre locations (<pi,<p2,4>3) — (60°, 180°, 300°). In general, dissymmetry of the current has 
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been found to be more pronounced when the absolute difference between the three gap widths (in this case, 9°) is 
increased. Similarly, loading the gaps with lumped impedances whose values are sequential in clockwise vs. anti- 
clockwise direction produces dissymmetric azimuthal currents, even at exceedingly low frequencies {ka <C 1) where 
distributed effects are small. Figure [5] shows an example where different lumped impedances Zk are realized by 
filling the gaps <50fc with different dielectric inserts of relative permittivities e/eo = 1> 10, and 100, arranged in either 
clockwise or anti-clockwise fashion. Similar dissymmetric current distributions are found at higher frequencies and 
for oblique incidence. 



ka=1.5, (<!>., <|>,,<L) = (60,180,300) deg 
4, x10- 3 , 




100 200 300 

<|> (deg) 



Figure 4. Real (solid and dashed curves) and imaginary (dash-dotted and dotted curves) parts of 
the current I(cj)) induced in a passive three-split loop with ka = 1.5, a — 20.5 mm, ro = 0.82 
mm. Calculation is based on the first twenty coefficients l/a n (n = 0, . . . , 19) for the closed loop 
and the coefficients /„ for illumination by an obliquely incident plane wave polarized along the y- 
direction, i.e., (6, <f>o,tp) = (30°,0°,0°), for a planar chiral split ring with (0i, 02,03) = (60°, 180°, 300°), 
(501,502,503) = (1°,10°,19°) and for its planar enantiomer with (01,02,0s) = (60°, 180°, 300°), 
(501,502,50a) = (19°,10°,1°). 



3. Multipole characterization 

3.1. General considerations 

The Fourier expansion method as applied to mirror multi-split loops allows for an investigation of the effect of 
planar chirality on the current distribution. The induced current is fundamental to the scattering properties and 
-when formed into arrays or gratings- to the reflection and transmission characteristics. This allows for microscopic 
multipole modelling of a single loop and calculation of macrocopic constitutive parameters for the effective medium 
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Figure 5. Real and imaginary parts of the current 7(0) induced in a passive geometrically symmetric 
three-split loop with ka = 0.05, a = 20.5 mm, r = 0.82 mm, (01,02,03) = (0°, 120°, 240°), 6 (pi = &<t>2 = 
(503 = 3°. Calculation is based on the first twenty coefficients l/a n (n = 0, . . . , 19) for the closed loop and 
the coefficients /„ for illumination by a normally incident plane wave polarized along the y-direction, i.e., 
{0,cj>o,ip) — (0 o ,0°,0°) for anti-clockwise-sequenced gap permittivities (e r i, e r 2, e r 3) = (1,10,100) and for 
the planar enantiomer (clockwise-sequenced) arrangement with (e r i, e r 2, £13) = (1, 100, 10). 



formed by a collection of such loops. The latter can be carried out with the aid of an appropriate mixing rule 
for composite materials, e.g., using mixing rules based on Maxwell Garnett, Bruggeman, or differential equation 
formalisms [30 . For an example of numerical calculation of the effective permittivity, permeability, and chirality 
dyadics of a homogenized metamaterial consisting of planar chiral spiral inclusions using the Maxwell Garnett 
formalism, cf. figures 8 and 9 of [2D]. The higher-order multipolarizabilities of each loop can be used to define 
higher-order macroscopic permittivity, permeability and magneto-electrical coupling polyadics [25 , with which the 
macroscopic effects of locality and spatial dispersion can be characterized in an explicit manner. 

3.2. Multi-polarizabilities 

In general, electric and magnetic fields at r' radiated or scattered by a particle of volume V with current density 
J(r) of harmonic time dependence exp(jwi) can be related to its multipole moments by a Taylor expansion [31j-|33] 
of the vector potential of J_(r), i.e., 




(39) 



For a thin- wire circular loop of length L carrying a local current £[r(r, </>, z)\ 



J(r) = iy>)5(r - a)8{z)lM). 



(40) 
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The electric and magnetic dipole moments of this loop current are |33] 

£c ^ (jc^)- 1 J£ J(r)dV = (jc)- 1 J l(i)M (41) 



and 



P m "y|^fc)xJfc)dl/ = ^ J qWxIWM, (42) 
respectively, where the second equalities in (|4Tj) and (|42|l rest on the thin-wire assumption with line coordinate i. 



These expressions apply with E_ and H_ = B_/no as the source fields. 

Dipole moments induced by an externally incident wave can be linked to the source fields via dipolarizability 
and multi-polarizability polyadics, as outlined in | Appendix A| 21j. For the present purposes, we may write (|A.1|) 
(|A.2|I symbolically as [34] 

P c = eogjY.) • + VM^^CV) ■ (43) 
P m = VM^l mo (Y) • E™ + Mol mm (Y) • £ inc , (44) 

where the product operators 7r, . ( V ) incorporate effects of spatial dispersion to any order, through their dependence 
on V and powers thereof, interpreted as polyadic products. The dipolarizabilities and higher-order multi- 
polarizability polyadics can be extracted, for example, using a recursive multi-polarizability algorithm (RMA) 
based on the method of counterpropagating waves [21, 25, 35j. For its application to planar structures, the method 
involves some subtleties which are outlined in | Appendix B| 

For an active split ring, the dipole moments are not produced by external illumination with an incident field. 
Its multi-polarizabilities depend instead on the applied voltage and spatial gradients of the resulting electric field. 
For example, for the electric field E_ k = Vk/ {84>k)'L c j >k generated by Vk applied across the gap S(f>k at (p — <pk, the 
magnetic dipole moment resulting from the impressed azimuthal current defines a magneto-electric dipolarizability 
dyadic as 

E m = VMM mc '^- (45) 
3.3. Electromagnetic signature of planar chirality 

As emphasized in [23], chirality is a purely geometric notion. While it manifests itself through symmetries and 
dissymmetries of the azimuthal current distribution, this does not necessarily affect electromagnetic (EM) and 
magnetoelectric (ME) polarizabilities through pseudo-scalarity of their elements. Indeed, the integrated azimuthal 
current J L I((f>)d(f) (which is fundamental to calculation of all mulitpolc moments) is itself insensitive to planar 
dissymmetries. By subtracting the bidirectional excitation from the unidirectional one, we obtain an EM signature 
of planar chirality of a structure, for a given direction of incidence and polarization. For normal incidence, this 
signature can be defined as 



±ke,af3 2 \—ke,a/3 + ±-kl,a(i ) ~ 2 \— ki < a —kl,a{3 



(46) 
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where the upper and lower signs refer to a planar chiral structure "(+)" and its enantiomer "(— )", i.e., for a 
planar chiral structure "(+)" and its enantiomer "(— )" or, equivalently, for propagation in +fc 7 - or — fc 7 -direction, 
respectively (a, /3, 7 = x, y, z\ 7 ^ oc,0). Several other measures of planar chirality have been proposed in recent 
years, cf., e.g., [21]-[24] and references therein. 

3.4- Numerical example 

Figure [6] shows elements of the electric-electric and magneto-electric polarizability dyadics 7r cc and 7r mc for a 
three-split ring in which the locations of two splits are fixed at (fa, fa) = (0°,180°) whilst the third split at 
fa is being moved in positive direction around the loop, from 0° to 360°. All splits are given equal widths 
Sfa = Sfa = Sfa = 1°. Quasi-static excitation of the split ring is considered here (ka = 0.05 <C 1; 
A(l3 — [Mo/(47r)] J, J[r(^)]/||r' — r(^)||d£). The series expansion for the current is limited to its first twenty 
Fourier coefficient terms. 

From the Figure, it is observed that 

7tee,yx — i ]^ee,yy, 

T^mc, zy — i ]^rac,zx 7 (^'v 

implying that the left-hand (LCP) and right-hand (RCP) circular wave polarization states E x ± ]E y are 
eigenpolarizations of the scattered field. Secondly, with this choice of fa and fa, two loops that have their third 
split at fa or —fa are enantiomers because fa and fa are diametrically opposite. Therefore, for this choice of 
locations for the fixed gaps, 7r ee ^ and 7r mGi 2* are pseudo-scalars (i.e., their signs become inverted upon taking the 
planar mirror image), whereas 7r CCja » are true scalars. The latter follows from the fact that the two fixed gaps are 
both located on the ox-axis. Thus, 

J+) - - 

ir {+) - 

_(+) _ _ „(-) fAsl \ 

where n$ a p = nkt,ap(<i>i, <f>2, ^3) and n$ a p = irkl,ap{-4>\,-fa,-<fo) with, in this case, fa = - fa and fa = - fa. 
Equivalently, it follows that Tr CCtX * and 7r mCiZ * (but not 7r CCiJ( », because of the particular choice of fa and fa) are 
signatures for planar chirality in this case and that LCP and RCP are eigenpolarizations. Obviously, the situation 
can be symmetrized with respect to y by choosing fa — ±45° and fa — ±225° instead. 

When the two fixed splits are not diametrically opposite, the azimuthal dependence of 7Vki, a ,0 is no longer 
symmetric, as witnessed from Figure [7] for the case fa = 120° with all other parameters left unchanged. In this 
case, the Wkl,y* are no longer true scalars but also satisfy (|48|) . Comparison with Figure [6] further shows that the 
change in the irke,a,f3 is most pronounced when fa is close to fa. The Figures also suggest that a jump caused by 
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a current discontinuity occurs whenever one split eclipses another one, although it has been found that the size of 
this discontinuity becomes smaller when the 5(4>k) are made smaller (not shown). A more detailed investigation 
using a higher number of Fourier coefficients would be needed to analyze this effect in more detail. 

4. Reflection and transmission for gratings of planar chiral split rings 

Full-wave numerical simulation using periodic method of moments 

Large collections of multi-split loops give rise to radiated, scattered, reflected and/or transmitted waves that can 
be characterized with the aid of (|39[) — (|40p provided I(4>) in such collections is known. Planar arrays of closed loops 
and other elements have a long history of use as frequency selective surfaces (diffraction gratings) for application 
in reflector antennas [361 137j . For such arrays, it is appropriate to use a full- wave numerical method for calculating 
induced currents (magnitude and spatial distribution), in order to account rigorously for mutual coupling (EM 
interaction) between the elements [20j . This mutual coupling causes the amplitudes of the Fourier expansion 
coefficients u(4>) and y(cf>) to be different from those for a single isolated split ring. For periodic arrays, an efficient 
rigorous method of analysis is the periodic method of moments (PMM) [38) . Its implementation based on subdomain 
functions is used here, which has previously been detailed and validated for an array of two-split rings [39j and 
other types of wire elements [3D] • 

The subdomain PMM formulation allows for easy calculation of mutual coupling between non-corresponding 
wire segments in different cells of the periodic array (off-diagonal elements of the impedance matrix). However, 
it is in principle possible to set up the entire impedance matrix for curved wire segments as well. Although 
we are focusing here only on specular reflection or transmission characteristics for normal incidence (involving 
neither diffraction, surface waves, nor grating lobes) at wavelengths that are large compared to the elements' size 
and spacing, we use a relatively large number of spectral plane-wave components n x x n y in the simulations for 
calculating the array impedance elements, viz., n x = n y — 10. Previously, we have shown [39j that small values 
of n Xl n y (in particular, n x = n y = 1 as would be intuitively chosen in the present case) may produce highly 
inaccurate and unstable values for the mutual impedance elements and, hence, for the reflection or transmission 
characteristics. 

In randomized arrays of (planar chiral) inclusions, interaction effects are implicitly taken into account within 
the framework of effective medium theory using mixing rules [20] . as mentioned before. Alternatively, the interaction 
between particles can be explicitly accounted for through calculation of the modified polarizability tensors |41j . 
which do not require one to choose a mixing rule but instead calculates the perturbation of the microscopic current 
density from the particular configuration. The latter method is usable for regular and irregular arrays, at arbitrary 
wavelengths, but becomes cumbersome when the number of particles is large. 
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4-- 2. Numerical results 

Figures |5^a)(b) and[Hta)(b) show frequency characteristics of the in-phase and quadrature components of the cross- 
and co-polarized transmission coefficients, respectively, for a periodic array of triple-split loops. The split rings 
form a regular lattice with period 49.2 mm in both x- and y-directions. The loop radius and wire radius are 
a = 20.5 mm and ro = 0.82 mm, respectively. The gap locations and split widths are listed in the Figure captions. 
Each of the four arcs, discretized using piecewise straight angled wire segments, was assigned four internal nodes 
for calculating the current distribution. The bandgap can be shifted to any desired wavelength (within the classical 
field model) by appropriate scaling of the geometrical parameters. The characteristic effect of the planar chirality 
[5] on the phase of the coefficients, viz., an inversion of the signs of r xy , r yx , t xy and t yx while maintaining the 
signs of r xx , r yy , t xx , and t yy is clearly witnessed. 

Figures Etc) and [9fc) show the normalized intensity of the cross- and co-polarized transmission, |^t' ) | 2 and 
|i^| 2 , respectively. The cross-coupling, although small, is seen to be anti-symmetric with respect to handedness 
(ty$ = —tyx) and is sharply peaked near the resonance frequency. The figures confirm that the specific sense of the 
planar handedness is invisible to the transmitted cross-polar (non-diffracted) intensity, \t y t}\ 2 . Figure [9] shows that 
the resonance frequency of the transmission coefficient for y- (but not for x-) direction exhibits a high sensitivity to 
planar chirality, as an indirect consequence of the different lengths of wire segments, although not to the particular 
sense of handedness. The effect can be used for highly-sensitive detection of asymmetries or dissymmetries with 
respect to specific reflection planes. Similar results were found to apply to the cross-polarized reflection coefficient 
(ryx — —fyx) and associated reflected intensity Ir^l 2 . 

It should be noted that the choice of the reference azimuth cf> — for planar chiral rings carries physical 
significance because such rings are plano-chiral, rather than plano-semi-chiral [5] like e.g. four-legged swastikas or 
other types of centro-symmetric planar structures. Because of this lack of a point centre of symmetry, circular 
and linear birefringence effects are coupled. The latter is apparent by comparing Figure [9] with Figure 1101 
Parenthetically, note that the second resonance of t X x exhibits a very high Q-factor. In gratings of planar chiral 
split rings, the linear birefringence of planar chiral rings may be suppressed by randomly orienting them in the 
plane, giving rise to randomization of the values of <fio £ [0, 2tt], cf. [20]. While the PMM analysis technique 
requires strict periodicity only in one direction, as a minimum [381 Sec. 3.2], aperiodicity in the second dimension 
would itself introduce a preferential direction inside the plane (array effect) and hence would not eliminate the 
possibility of linear birefringence. Therefore, only a more cumbersome standard method of moment analysis [42] 
allows for analysis of arrays of randomly rotated planar chiral rings to resolve this problem. 

5. Conclusions 

Asymmetrically split rings serve as canonical examples of planar chiral structures, which allow for discrimination 
between circular polarization states of opposite handedness by a phase-coherent detector. In this paper, we 
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investigated the effect of planar handedness on the azimuthal current distribution and polarizability coefficients of 
a single multi-split ring, and the reflection and transmission properties of regular 2-D arrays of such rings. It was 
found that geometrical as well as EM constitutive properties of the rings and their gaps allow for dissymmetric 
characteristics by introducing clockwise vs. anti-clockwise progressions of gap widths or permittivities of gap inserts, 
respectively. An analysis method based on the Fourier expansion of the current was introduced that provides an 
efficient alternative to full-wave simulation of an active or passive multi-split ring. 

Increasing the number of gaps may be expected to give rise to more complicated frequency characteristics, 
because of the interplay of resonances associated with each of the wire segments whose number increases with 
increasing number of gaps. This requires an increasingly large number of Fourier coefficients of the wire current to 
be retained in order to maintain numerical accuracy. 

Acknowledgments 

The author wishes to thank the reviewers for their constructive comments and suggestions. 



Planar Chirality of Plasmonic Multi-Split Rings 



19 




Figure 6. Real part (blue solid curve) and imaginary parts (green dashed curve) of selected polarizability 
elements 7r COiQ/ 3 and 7r mCiZ( 3 (in units 10 -3 m 3 ) for a single three-split ring for external plane-wave quasi- 
static illumination (ka = 0.05) with (0i,0 2 ,0 3 ) = (0°, 180°, 0° . . . 360°) and Sfa = &4>2 = 8<j) S = 1°: (a) 
for the original structure with incidence in negative z-direction, and (b) for the enatiomeric structure (i.e., 
for incidence onto the original structure in positive z- direct ion ) . 
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Figure 7. Same as Figure [6] but for (a) <j) 2 = 120° and (b) <j) 2 = -120°. 
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Figure 8. (a) In-phase component, (b) quadrature component, and (c) squared magnitude of cross- 
polarized transmission coefficient t yx for a 2-D periodic array of triple-split planar chiral loops with 
(0i, 02,0s) = (0°, 120°, 240°) and (50i, 50 2 , 50 3 ) = (5°, 5°, 5°), (3°, 5°, 7°), or (3°, 7°, 5°). Corresponding 
values for are identical. 




123456 123456 123456 

Frequency (GHz) Frequency (GHz) Frequency (GHz) 



(a) (b) (c) 

Figure 9. (a) In-phase component, (b) quadrature component, and (c) squared magnitude of co- 
polarized transmission coefficient t yy for a 2-D periodic array of triple-split planar chiral loops with 
(01,02,0a) = (0°, 120°, 240°) and (501,502,503) = (5°, 5°, 5°), (3°, 5°, 7°), or (3°, 7°, 5°). Curves for 
enantiomeric forms are indistinguishable from each other. 
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Figure 10. (a) In-phase component, (b) quadrature component, and (c) squared magnitude of co- 
polarized transmission coefficient t xx for a 2-D periodic array of triple-split planar chiral loops with 
(4>i,<fa,<j> 3 ) = (0°, 120°, 240°) and (6<t>i, 5<j> 2 , <%) = (5°, 5°, 5°), (3°, 5°, 7°), or (3°, 7°, 5°). Curves for 
enantiomeric forms are indistinguishable from each other. 
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The multi-polarizabilities of a point-polarizable particle can be expressed by means of a Taylor series expansion 
for its electric and magnetic dipolc moments, as defined by ([4"T]) and (|42[) . with respect to spatial derivatives of the 
excitation (source) fields as 125, |33] [43] 
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in which nil, single, double, triple, etc., underlined quantities represent scalars, vectors, dyadics, tryadics, etc., i.e., 
tensors of rank zero, one, two, three, etc., respectively. We limit the further analysis to first-order terms in V. The 
auxiliary tryadics q' are related to the actual quadrupolarizabilities q via [33] 

q' =±q T ^ 1 (A.3) 

=ki ^lk 

where the upper and lower signs correspond to k = £ and k ^ £, respectively, and where the superscript T^ 



denotes backward monadic transposition, i.e., 
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moments q and q , here defined by [32] Ch. 2], [33 
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the quadrupolarizabilities q are obtained from 
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+ 0(V 2 ), (A.7) 
when I(£) is known, where o'W = ±oW . By extension, the octupolarizabilities o follow from a similar expansion. 
For a density of JV particles per unit volume, macroscopic multipole moments and multi-polarizability 



polyadics, denoted by corresponding upper-case symbols, can be defined through P_ c = Np^ = £oP. cc ' He + 
\/lM^oE. cnl ' ii'i Q = Nq('\ etc., in which the original external incident plane- wave fields P mc and H_ mc can 
be replaced by the local internal fields and IP to account for mutual interaction, with which they coincide 
in the dilute limit. The microscopic multipole expansion of p and p m can be extended to a macroscopic multi- 
polarizability expansion for and P m and, hence, for e and /i [25] . This allows for identifying, explicitly, the 
contribution of spatial dispersion in the permittivity and permeability dyadics. 

Appendix B. Method of counterpropagating waves for extracting multi-polarizabilities of planar 
structures 

The method of counterpropagating waves [35] can be used for the extraction of dispersion characteristics of all 36 
elements of the four 3x3 electro-electric, magneto-magnetic, electro-magnetic and magneto-electric constitutive 
dyadics. The method has been applied to previously to 3-D chiral inclusions (helix [25], composite twisted hooks 
[34], etc.). For analysis of 2-D planar structures, its application requires due care: waves that are incident from 
opposite sides experience enantiomeric forms of such structures. Therefore, the local negative superposition of 
waves and dipole moments is, in general, not valid in this case. However, for normal incidence, the magnetic field is 
parallel to the plane of a loop and hence will not induce currents. A similar situation occurs for grazing incidence 
when the electric field is perpendicular to this plane. For example, the excitation (E x , H y , +k z ) characterized by 
(9 = 1 80°, 0o = —90°, V> = 0°) yields two of the electric-electric dyadic elements: from p C:X = eon cc _ xx E x (+k z ) and 
(|3T]). it follows that 



]uje E x {+k z ) 

; o 2 V(0)sin(0)d0 (]rn 



juje E x (+k z ) 

ee,vx juje E x (+k z ) 

= \T 7 <» cos(0)d0 

ju€ E x (+k z ) ' 1 ■ ' 

where /(</>) = J(0)1j,(0) is the azimuthal current flowing in the wire filament at 0, whereas for propagation in 

y-direction, using the excitations (E x , —H z , +k y ) and (E x , H z , —k y ) characterized by (9 = 90°, 4>o = 90°, ip = 180°) 

and {9 = 90°, 0o = -90°, V = 0°), respectively, 

r +) = 1 ( fqiW-U fUWM-W 

cc ' xA V) 2 \juje E x (+k y ) + ^e Q E x {-k y ), 
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(B.3) 

ccyAVy) 2 \^E x {+k y ) + )Lue E x (-k y )J ' 

(B.4) 



Likewise, from p m>z = y/noe n mCiZX E x (+k z ) and (g2|), or from p c<x = y/fXoeoTT em , xy Hy(+k z ) and |gT]), we obtain 
(+) (y? , [j^J L g(e)xI(£)d£-l z 



"'emliv(^z) 



e 2 E x (+k z ) 

e 2^(+fc z ) 1 ^ 



jw^oeo H y (+k z ) juje E x (+k z ) 

For the reversely propagating excitation (E x , — H y , — k z ), corresponding to (6 = 0°,cf>o = — 90°, tp = 0°), the 
obtained p cc _ XXl p cc , V x, and p mc ,zx are again those for the enantiomeric clement, which we denote by the superscript 
"(-)". Thus, instead of (|B~T]) , (IB~2|) . and JR5J, we now have 



n (-) (_ Vz ) = ^W d< ^ = jT^ww 

cc,xx Z )LoeoE x (-k z ) }LueoE x (-k z ) 



(-) r-V 1 = Jl^^ = / o 2 ^(0)cos(0)d^ 
^"'"^ " J jcje E x (-k z ) jue E x (-k z ) 



(B.7) 



"mc.zil v 2 



(B.8) 

^ J L #x/(<)dM 2 

e 2 
^ a/ o 2 "/(0)d0 

e 2E x {-k z ) ' 1 ' j 

Note that the superposition of bidirectional excitations, [i? a; (+fc !y ) + E X (— k y )]/2, produces a spatially harmonic 
standing wave in the direction of the countcrpropagation having sinusoidally varying amplitude with maximum 
located at y = 0, whereas (E x , H y , +k z ) represents a travelling plane wave whose amplitude is constant along the 
y-direction. For quasi-static excitation (ka <C 1), both excitations yield identical dipolarizabilities. For h ^ 1, 
however, the results can be substantially different. Countcrpropagation has the advantage that any contribution by 
H z , whether significant or residual, is completely cancelled, but has the disadvantage that the effect of the standing 
wave (i.e., contribution of spatial dispersion via the gradient of E x in y-direction and the quadrupolarizability 
elements q' ee ^ yx ) must be compensated for. Specifically, 

- (41,L - q -^v y ) Ex (b.io) 

pityEy = UV, ay ~ %V.) E y (B.ll) 
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,(+) 



tt {+) --(a (+)l V +o (+) ' V ^ 



E z 



(B.12) 



etc., where fc, £ = e, m; a = x, y, z; and with the spatial derivatives given by 



1 



1 



~VyE x EE - [VyE x ( + ky) + V v E x {-k y )] 

1 d 

= 2^[° x P( _ j fc «y) + cx P(i k yy)]Ex{±k y ) 
= - k y sm(k y y)E x (±k y ) 
etc., for —a < y < a. The respective dyadics Vfi and V if are given in full by 

f —k x sm(k x x)E y (±k x ) 



(B.13) 
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(B.14) 
(B.15) 



(i -fc-r sxn(k x x)H y (±k x ) -k x sin(k x x)H z (±k x ) ^ 

-k y sm(k y y)H x (±ky) -k y sm(k y y)H z (±k y ) 

i sin(/c 2 z)i7 x (±fc z ) — fc z sin(fc z z)i/j / (±/c 2 ) y 

Through cyclic permutation, piJLy, Pcc^j/, and Pme,zy are obtained from excitation by (E y ,—H x ,+k z ) given by 
(9 = 180°, <po = O ,^ = 0°). For the enantiomeric form, pi e } xy , pie}yy, and piae,zy follow from excitation by 
(E v , H x ,~k z ), i.e., (9 = 0°, O = 0°, if, = 0°). 

For the remaining electro-magnetic and magneto-magnetic dyadic elements, excitation by (E x , H z , —k y ) can be 
superimposed with a matching excitation (—E x ,H z ,+k y ), which correspond to (9 = 90°, (j>o — — 90°, -0 = 0°) and 
(9 = 90°, 0o = 90°, ifj — 0°), respectively. Similarly, excitations produced by (E y , H Zl +k x ) and (—E y ,H z ,—k x ), 
corresponding to (0 = 90°, O = 0°,V> = 0°) and (9 = -90°, <f> Q = 180°, = 180°), can be superimposed. Thus, 
from 2p CiX = yfjm^p cm , xz [H z (-k y ) - H z (+k y )], 2p miZ = HoPmm,zz[H z (-k y ) - H z (+k y )], and (gT|), we obtain 
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(B.16) 
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in which H z (±k x ) and H z (±k y ) refer to the magnetic field associated with waves propagating in diametrically 
opposite directions. If the field components H a are obtained from the same set of counterpropagating waves as 
for E ai i.e., now through their negative superposition, then again a standing wave results. This standing wave 
is now in quadrature with respect to the one for the electric-field excitation and extends along the direction of 
propagation, e.g., 

1 

2 



- [V z H y (+k z ) - V z H y (-k z )} 



1 d 

= r^- [cxp(jfc z z) - exp(-jfc z z)] H y 
z oz 

= jk z cos(k z z)H y , (B.19) 

which, for the planar structure, is to be evaluated at z = 0. Note that a 2-D planar element (irrespective of whether 
it is planar chiral, symmetric, or asymmetric) exhibits nonvanishing spatial dispersion in this plane because of the 
nonzero spatial extent of the element in its plane. Specifically, for grazing incidence, E x (resp. E y ) traverses a 
finite length of the element during propagation in ±y- (resp. ±x-) direction and thus contributes to p m ,z- This 
does not occur for these polarizations when propagation is in indirection, because of the zero thickness of any 
2-D element in this direction. Similarly, H x and H y affect the induced overall p e for grazing incidence (±k y or 
±fc x , respectively). 

The elements of the quadrupolarizability tryadics are obtained by extension: from (|A.5|) and with the same 
plane-wave excitations as for the dipolarizabilities, we obtain, up to first order in V, 

(+ ) _ kW)m+mw)w--uh m on , 

a (+) 

3 E x (+k z ) 

(B.21) 

etc. These allow for computation of first-order corrections of the dipolarizabilities for spatial dispersion effects. 



